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ABSTRACT 

We use two-dimensional kinematic maps of simulated binary disc mergers to investi- 
gate the A/j-parameter, which is a luminosity weighted measure of projected angular 
momentum per unit mass. This parameter was introduced to subdivide the SAURON 
sample of early type galaxies in so called fast Xn > 0.1 and slow rotators Xji < 0.1. 
Tests on merger remnants reveal that Xn is a robust indicator of the true angular 
momentum content in elliptical galaxies. We find the same range of Xn values in our 
merger remnants as in the SAURON galaxies. The merger mass ratio is decisive in 
creating a slow or a fast rotator in a single binary merger, the former being created 
mostly in an equal mass merger. Slow rotators have a which does not vary with 
projection. The confusion rate with face-on fast rotators is very small. Merger with 
low gas fractions form slow rotators with smaller ellipticities and are in much better 
agreement with the SAURON slow rotators. Remergers of merger remnants are slow 
rotators but tend to have too high ellipticities. Fast rotators maintain the angular 
momentum content from the progenitor disc galaxy if merger mass ratio is high. Some 
SAURON galaxies have values of A^ as high as our progenitor disc galaxies. 

Key words: methods: analytical - methods: N-body simulations - galaxies: elliptical 
and lenticular, cD - galaxies: formation - galaxies: evolution - galaxies: fundamental 
parameters - galaxies: kinematics and dynamics 



1 INTRODUCTION 

The connection between projected and the intrinsic prop- 
erties of elliptical galaxies is a notorious problem. The 
Hubble classification of early-type galaxies is based on 
their elliptical shape, which was the natural observable 
distinctive feature. It was already pointed out by Hubble 
that this sequence is only a sequence of apparent shape 
depending on the inclination under which we observe a 
galaxy. When kinematic observations of ellipticals became 
available, they showed that some of them rotate t o o slowly 
to be shaped by rotation alone ()lllingworthlll977h .lB mnevi 
() 19781 ) proposed that these systems are probably shaped to 
a large extend by anisotropic pressure. This is best seen 
when the apparent flattening is plotted versus the balance 
between ordered and unordered motion {v/a,e) using the 



relation of these observables for an isotropic rotator as a 
reference point. Observations also showed that the am ount 
of rotation varies with luminos ity l|Davies et al.|[l983 ) and 
isophotal shape (|Bendeij[l988l ). Isophotal shapes as mea- 
sured by the 04 parameter can either appear boxy (04 < 0) 
or discy (04 > 0). Boxy and discy galaxies seem to form 
a true dichotomy as the isophotal shape correlates with 
on first sight unconn ected properties like X-ray luminosity 
Bender et al. 19891 ) and central surface density slope 



Ferrarese et al .' 'l994':'Laucr et al.' 'l995l lFaber et al.l [1993). 



This led iKormen dv & Bender (1996) to propose a revision 
of the Hubble sequence using the isophotal shape to cast el- 
liptical galaxies into boxy and discy ellipticals, which might 
be a classification more closely con nected to their formati on 
history than apparent ellipticities (|Naab fc Burker3l2003l ). 
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Kinematic properties are normally measured from long- 
slit data at the major axis (v) and a central aperture (a). 
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But these are very special positions, which might or might 
not be representative for the overall structure of the galaxy. 
It is now possible to extract the line-of-sight velocity dis- 
tributions (LOSVDs) from the full two-dimensional field 
of a galaxy with modern integral-field instruments. The 
SAURON survey pioneere d the first comprehen sive sample 
of nearby E/SO galaxies (|de Zeeuw et al.l I2OO2I ) with two- 
dimensional kinematic data. Some of these galaxies exhibit 
very complex velocity or velocity dispersion fields which are 
not easily captured b y the usual v/a,e parameterization 
l|Emsellem et aLlbOO^ ). It can happen that two galaxies with 
similar long-slit measurements of ellipticity and v/a have 
morphologically different kinematic fields, e.g. one galaxy 
has a kinematically de coupled compon e nt an d the other 
a regular velocity field. lEmsellem et al.l (|2007l ). henceforth 
EM07, therefore introduced a luminosity weighted measure 
of the specific line-of-sight angular momentum 



< R\V\ > 
< RVV^ + a2 >' 



(1) 



where the pointed brackets are luminosity averages. is 
designed to account for radial changes of angular momen- 
tum content in kinematic fields. T his is dif f erent to the mea- 
sure < V > / < a > proposed bv iBinnevI (|2005| ) to ex ploit 



integral field kinematic data. However, BinnevI l|2005h also 
stressed that with integral field data available, the connec- 
tion between internal and projected structure of a galaxy is 
much less ambiguous. 

The E/SO galaxies are cast into the category of slow ro- 
tator or fast rotator according to their measured A_r (EM07): 
slow rotators ha v e a X r < 0.1 and fast rotators Xr > 0.1. 
ICappellari et al.1 (|2007l ) showed that the slow and fast ro- 
tating galaxies of the SAURON sample populate distinct 
regions in the (n/cr, e) diagram, even if their inclination is 
taken into account, which means that those galaxies are also 
intrinsically fast or slowly rotating. Slow rotators tend to be 
more massive, have lower ellipticities and exhibit faster ro- 
tating, kpc-sized, old kinematically decoupled components 
(KDCs), while fast rotators have lower luminosities and can 
also have central KDCs, which are howe ver small and of- 
ten c onsists of young stellar populations (jMcDermid et all 
l2006h . Slow rotators also have strong kinematic twists and 
kinematic misalignments, while photometry and velocity 
contours are almost always aligned in fast rotators. In ad- 
dition, a correlation between anisotropy and ellipticity has 
been detected for slow and fast rotators with slow rota- 
tors being more isotropic and fast rotators b eing anisotropic 
ICappellari et al.ll2007l : iBurkert et ahlbOOTh . 

In this paper we want to connect the findings of the 
SAURON ga laxies with their fo r matio n history. Ever since 
the work of iToomre fc Toomrd (|l972l ') the merger of two 
late-type galaxies has been so far the most appealing and 
the best tested theoretical model for the formation of an 



elliptical Ealaxv (BarnesI 19921: Barnes & Hernauist 199(51: 


Naab & BurkertI 


2003 


; Cox et al.l 20061: Robertson et al.1 


20061: iNaab et al.1 


2006 


and references therein). 



The full two-dimensional kinematic field of merger 
remnants have sca rcely been studi e d. Fol lo wing the pi- 
oneeri ng work of iBendo fc BarnesI (|2000l ). Ijesseit et al.l 
l|2007t l. henceforth J07, analyzed kinematic maps of N-body 
merger remnant s with kinemetric methods, introduced by 
iKrainovic et al.l (|2006l ). They quantified the presence of 



counter-rotating cores, the fiattening of isovelocity contours, 
counter-rotating discs and other peculiar features which are 
also observed in the SAURON galaxies. In general, maps 
of slow rotators resembled more equal-mass mergers while 
fast rotators show kinematical fields similar to unequal-mass 
mergers. However, with the introduction of Xr, we can quan- 
tify the line-of-sight angular momentum content of N-body 
merger remnants in the same way as has been done with 
observed galaxies. Ideally we can connect merger remnants 
which appear as slow rotator or fast rotator with a certain 
formation mechanism. Beyond that we know the true angu- 
lar momentum content of a simulated N-body merger rem- 
nant and can therefore determine whether the remnant is 
indeed intrinsically a slow (fast) rotator or not. 

The aim of this work is two-fold: we want to study how 
much can be learned from the intrinsic structure of a galaxy 
by calculating its line-of-sight angular momentum parame- 
ter Xr. In addition, we want to clarify whether remnants 
of binary mergers follow the slow-fast rotator division of 
the red sequence as found in the SAURON sample. The pa- 
per is organized as follows: we describe the simulations in 
Sec[2l give details about how we calculate the maps and de- 
termine A_R (Sec l2.2|) and examine the connection between 
line-of-sight angular momentum and intrinsic angular mo- 
mentum (Sec[3}, show the projectional behaviour of Xr for 
merger remnants with varying orbital content (Sec|4]), clas- 
sify the maps derived from the merger remnants as fast or 
slow rotators and investigate the influence of various forma- 
tion mechanisms (Sec[5} and finally discuss the results in 
Seed 



2 SIMULATIONS 

2.1 The Merger Models 

The model parameters of the sample of 1:1 and 3:1 disc 
mergers we are studying in this paper is identical to the 
sample presented in NJB2006, except that we include star 
formation and stellar feedback in the modeling of the dissi- 
pative component. Here we give a brief review of the setup. 

The disc galaxies were constructed i n dynamica l equi - 
librium using the method described by iHernguisO (|l993l ) 
with the following system of units: gravitational constant 
G=l, exponential scale length of the larger progenitor disc 
hd = 1 (the scale height was = 0.2) and mass of the 
larger disc Md = 1. The discs were exponential with an ad- 
ditional spherical, non-rotating bulge wi th ma ss Mb = 1/3, 
a Hernquist density profile (Hcrnguistl Il990l ) and a scale 
length r(, = 0.2/id, and a pseudo-isothermal halo with a 
mass Md = 5.8, cut-off radius Tc ~ lOhd and core radius 
"/ — Ih. The parameters for the individual components were 
the same as for the coUisionless mergers presented in NB03 
and NJB2006. We replaced 10% of the stellar disc by gas 
with the same scale length and an initial scale height of 

hz,gas — O.l/lj. 

Whereas the NB2003 and NJB2006 simulations 
had been performed with the Nbody- SPH code VINE 
llNelson et al.ll2008l : IWetzstein et al.ll2008|') we u sed the par- 
allel TreeSPH-c ode GADGET-2 (ISpringellfioOsh to follow ra- 
diative cooling (|Katz et al.lll99d ) for a primordial mixture of 
hydrogen and helium together with a spatially uniform time - 
independent local U V background (|Haardt fc Madaulll996l ) . 
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In addition, we include star formation and the associated 
supernova feedback followins the sub-resolution multiphase 
model developed bv lSprinee l & Ifcrnauist (2003') using the 
same implementation as inljohanssori ct al. (200^) but with- 
out feedback from black holes. 

We followed mergers of discs with mass ratios of 1:1 and 
3:1. The equal-mass mergers were calculated using in total 
440000 particles with each galaxy consisting of 20000 bulge 
particles, 60000 stellar disc particles, 20000 SPH particles 
representing the gas component in the disc, and 120000 halo 
particles. Twice as many halo particles than disc particles 
were used to red uce heating and instability effects in the 
disc components (|Naab et aLlll999l ) by encounters between 
halo and disc particles. For 3:1 mergers the parameters of the 
more massive galaxy were as described above. The low-mass 
companion contained a third the mass and the number of 
particles in each component, with a disc scale- length (stars 
and gas) of /i = -^1/3. 

The gravitational forces were softened with a Spline ker- 
nel of /igrav = 0.05. The minimal size of the Spline kernel 
used for computing the SPH properties, /isph, was fixed to 
the same value. Implicitly, this procedure suppressed gas col- 
lapse on scales smaller t han the softening scale and prevents 
numerical instabilities l|Bate fc Burkertj[l997l ). The initial 
discs were run in isolation for two dynamical times to al- 
low the systems to finally settle into an equilibrium state. 
The galaxies merged on parabolic orbits with a pericenter 
distance of 2 disc scale lengths. The initial geometries are 
identical to NB2003, Table 1. and NJB2006. The merger 
remnants were allowed to settle into dynamical equilibrium 
for approximately 30 dynamical time-scales after the merger 
was complete. Then their equilibrium state was analyzed. All 
simulation have been run on a SGI Altix at the University 
Observatory in Munich. 



profile using Gauss-Hermite basis functions (|Gerhardlll993l : 
Ivan der Marel fc Franx|[l993 ) . The kinematic parameters of 
each profile (o"flt, Uflt, h^, hi) were then determined by least 
squares fitting. 

\r can be calculated from integrated field data via the for- 
mula given in EM07 



(2) 



where Fi is the fiux,i?i the projected radius, Vi the line- 
of-sight velocity and ai the line-of-sight velocity dispersion of 
each grid cell. We determined these properties as explained 
in the previous section and can therefore calculate \r di- 
rectly from our 2D data. As galaxies are projected randomly 
on the sky, we project each remnant 200 times randomly on 
the sphere of viewing directions. This gives us a sample of 
20000 2-dimensional maps. SAURON has a finite field-of- 
view which extends to one effective radius, which we need 
to take into account as \r can vary strongly with R. We 
determine for each projection the effective radius and sum 
only over the grid cells inside one effective radius. Thus we 
ensure a fair comparison to the SAURON data. 



3 THE INTRINSIC SPECIFIC ANGULAR 
MOMENTUM AND \r 

3.1 Intrinsic Angular Momentum 

The total intrinsic specific angular momentum Jintr of a 
galaxy (or N-body merger remnant) can be calculated from 
the discrete stellar (particle) distribution. The components 
of J in Cartesian coordinates are 



(3) 



2.2 Determination of \r 

The two-dimensional maps examined in this article, are iden- 
tically determined as laid out in paper J07. We repeat here 
briefly the most important steps: For the 2D analysis we 
binned particles within the central 6 length units (21 kpc) 
on a grid of 48 x 48 cells. This correspond s typically to 2-3 
effective radii depending on projection (see lNaab fc Truiillol 
I2OO 6 for the exact determination of re). To include seeing 
effects we created for every luminous particle 10 x 10 pseudo- 
particles with identical velocities on a regular grid with a 
total size of 0.125 unit lengths (0.44 kpc) centered on the 
original particle position. The mass of the original particle 
was then distributed to the pseudo-particles weighted by a 
Gaussian with a standard deviation of 0.1625 unit lengths 
(0.56 kpc). Thereafter the pseudo-particles were binned on 
a 48 X 48 grid. For the kinematic analysis we binned (mass 
weighted) all pseudo-particles falling within each grid cell 
in velocity along the line-of-sight. The width of the veloc- 
ity bins was set to a value of 0.1 (26.2 km/s) for line-of- 
sight velocities vias in the range —4 ^ v\oa ^ 4 (corre- 
sponding to ±1048 km/s). This resulted in 80 velocity bins 
over the whole velocity interval. Using the binned velocity 
data we constructed line-of-sight velocity profiles (LOSVD) 
for each bin of the 2D grid. Subsequently we parameter- 
ized deviations from the Gaussian shape of the velocity 



(4) 



(5) 



and Jintr is the norm of that vector. For simplicity we implied 
here that all particles (stars) have the same mass and do not 
write the sums over masses explicitly. The angular momen- 
tum vector can projected to the plane of the sky, defining 
the projected intrinsic angular momentum vector Jp with 
components 



E 



E 



(6) 



(7) 
(8) 



where we define the y-direction arbitrarily as perpendicular 
to the plane of the sky. As our line-of-sight is along the 
y-direction we can not observe transversal motions, i.e. in 
X and z-direction. We can now define an apparent angular 
momentum vector Japp with the components 



Ja; — 



E- 



(9) 
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(10) 

(11) 



It has been show by iFranxl l|l988l ) that for triaxial galax- 
ies the apparent angular momentum amounts to half the 
intrinsic angular momentum, which we rewrite as 



Jp — KjJb 



(12) 



where kj is typically 2. An additional complication is that 
the angular momentum of a galaxy can be due both to 
streaming m otions and to figure rotation. Again we follow 
iFranxl (|l988l ) 



Japp — O.SJp 



i-projL^proj , 



(13) 



angular velocity of the projected density distribution, re- 
spectively. Indeed figure rotation is present in about 15% of 
the merger remnants, but we do not try here to disentangle 
the relative contributions of figure rotation and streaming 
motion. Although such a decomposition would be possible 
in principle for the N-body remnants such an analysis would 
be beyond the scope of this paper. 
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Figure 1. The apparent angular momentum calculated from par- 
ticle data (y-axis) as compared to the < R|V| > calculated from 
the velocity maps (x-axis). The agreement is very good for all 
projections and all merger remnants. 



3.2 Angular Momentum in Observations 

We cannot directly access the apparent angular momentum 
of a galaxy observationally. The closest proxy to the ap- 
parent angular momentum which we can calculate from the 
two-dimensional velocity field is < R|V| >, where R is the 
projected radius, V the line-of-sight-velocity and brackets 
denote luminosity averages (see also Appendix A of EM07 
for a discussion) . It is straightforward to calculate Japp from 
the N-body simulation and < R|V| > from mock observa- 
tions of the same remnants and same aperture. In Fig[T]we 
show that agreement is very good and verifies that the map- 
ping of the velocity fields gives a close account of the intrinsic 
structure. 

The inclination of a galaxy is in general not known and 
dynamical modeling o r other observatio nal clues are needed 
to infer its projection (|Cappellarill2008l ) . We can choose the 
projection of an N-body remnant at will and can calculate 
Japp and Jp for each projection. The ratio of these two values 
should be equal to two, if and only if we have the information 
over the whole field of the galaxy. This is in general not the 
case, as the field-of-view is limited by the instrument we 
are using. We are calculating therefore Japp and Jp for two 
different apertures, one which encloses the whole remnant 
and one which extends only to one efi'ective radius. This is 
illustrated in the right plot of Fig[2] where we highlighted 
the aperture at one effective radius. 

It is instructive to examine first a test model, for which 
we chose the progenitor disc galaxy, which is by construction 
axisymmetric and fast rotating. For almost all projections 
we find that kj (Fig[2] left) is very close to the theoretical 
value of 2, if calculated over the whole remnant. If we deter- 
mine Jp and Japp inside one Re we tend to get a spread of 
higher values for kj. This is not too surprising, as we catch 
the particles only on part of their trajectory (close to peri- 
center) and line-of-sight velocities are not representative of 
their total angular momentum content. This will depend on 



the detailed orbit structure of the remnant, i.e. how much 
angular momentum sits at large radii. In the same Fig[2]we 
illustrate the spread for the whole remnant sample, which is 
somewhat larger than for our axisymmetric test models, but 
with the same tendency to underestimate Jp. Again using 
the information on the full remnants shows that the appar- 
ent angular momentum and the projected intrinsic angular 
momentum are on average connected by a factor of two. 

The quantity Xr, as defined by EM07, has been pro- 
posed as a proxy for the intrinsic specific angular momentum 
of a galaxy within the radial range over which A_r is deter- 
mined. In the last test we want to see how closely Xr and 
the angular momentum content are connected for a random 
projection. We discussed before that < R|V| >, as a proxy 
for Japp, is closely connected to the projected intrinsic angu- 
lar momentum. Xr is in contrast to < R|V| > a normalized 
quantity. Therefore we scale both quantities to their max- 
imal values in order to compare the effect of inclination. 
We use again the progenitor disc galaxy as a test model 
for which ellipticity maps inclination accurately. We see in 
FigOthat Xr decreases much less rapidly with inclination 
than < RjVI > does, which is also shown by the histograms 
of the deviations from their maximal values (in this case 
scaled to one). The reason is that < v > and < a > de- 
crease simultaneously and their ratio stays roughly constant 
until we reach inclinations close to face-on. The true an- 
gular momentum content is of course most closely related 
to the maximum value of Xr, or < R|V| >, which in many 
cases, at least for fast rotators is a projection close to edge- 
on. In summary, even if there is a close connection between 
< R|V| > and the projected intrinsic angular momentum 
Jp , the projectional behaviour of Xr is more robust and will 
stay closer to its maximal value for a higher fraction of all 
possible viewing angles. 
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Figure 2. Left: Distribution of the ratio kj = Jp/Japp for the axisymmetric progenitor for 200 projections. The ratio shows the 
theoretically expected value of 2, if we calculate the angular momentum terms for the whole remnant, while a larger spread of kj is 
observed if we constrain the field-of-view to one Re- Middle: The same plot but this time for all coUisionless remnants. Right: illustration 
of the size of the aperture at one effective radius (white rhombus). 




Figure 3. Left: Xn, < R[V| > and Japp for random projections of the progenitor galaxy. All quantities are normalized to their maximum 
value. Right: Distribution of absolute deviation of the same quantities from their (normalized) maximum values. 



4 INTRINSIC SHAPE VARIATIONS: 
EXAMPLE REMNANTS 

Before we examine the properties of the merger remnant 
sample as a whole it is instructive to assess the variety of 
remnant properties present in the sample. The remnants in 
our sample have very different orbital makeups and intrinsic 
shapes. They can be dominated by major axis tubes, minor 
axis tubes and box orbi ts, and have prolat e, oblate or triax- 
ial shapes respectively jjesseit et al.l l2005h. Their projected 
properties must correlate with their intrinsic shape and it 
is interesting to study how Xr is varying with inclination. 
In Thomas et al. 2007 we selected merger remnants at the 
extreme ends of the shape distribution and we want to study 
the projected 2-dimensional properties of three of them to 
elucidate systematic kinematical trends with intrinsic shape. 

Our test remnants are t he OBLATE, PROL ATE 
and TRIAXIAL remnants from iThomas et al] l|2007l ). We 
are plotting the ellipticity, A_r and the isophotal shape- 
parameter a4 in Figure |4] to illustrate the variation with 
projection. The OBLATE remnant (red symbols) is the re- 



sult of a 4:1 merger, which is dominated by minor axis tubes. 
It is classified as a fast rotator from almost any viewing an- 
gle, although a dependence on inclination is visible. The 
TRIAXIAL remnant is a 1:1 merger remnant which is dom- 
inated by box and boxlet orbits. Interestingly the remnant 
has a Xh consistent with from all viewing angles. The 
PROLATE remnant (black symbols) has the most complex 
dependency on inclination angle. This is so because it con- 
tains a significant amount of both minor and major axis 
tube orbits. Therefore we almost always see some rotation, 
originating from different kinds of orbit classes (for the de- 
tailed de scription of the orbi t classification in N-body rem- 
nants see iJesseit et all (|2005h . However, as this is also a 1:1 
merger remnant it does not rotate very fast and reaches a 
maximum Xr of 0.18. 

The statistical distribution of Xr and e are better seen 
in the right column of FigH The Xr of the TRIAXIAL rem- 
nant has a delta-function like distribution at Xr = 0, that 
means it is a slow rotator under any viewing angle, while the 
OBLATE remnant is a fast rotator for almost all projections 
except for three inclinations close to face-on. This results in 
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a probability of « 1% to misclassify the OBLATE remnant 
as a slow rotator. The PROLATE remnant is one of the few 
remnants which could be classified as a fast or a slow ro- 
tator with equal likelihood, because both tube orbit classes 
seem to carry an equal amount of rotation and move along 
the line-of-sight for different viewing angles. The PROLATE 
remnant is rather round for almost all projections, while the 
TRIAXIAL and the OBLATE remnant are quite flattened, 
albeit the flattening is caused by different orbit classes, i.e. 
box orbits in the case of the TRIAXIAL remnant. Flattened 
slow rotators are observed in the SAURON sample only in 
exceptional cases (see also Sec[5]) 

4.1 Kinematical Misalignment 

Another discriminatory property of slow and fast rotators is 
the kinematical misalignment. The misalignment angle ^ is 
defined as the angle difference between the photometric and 
the kinematic position angle 

* = IPAphot -PAkinl, (14) 

where PAkin is defined as the angle 0° < PA kin < 180° along 
which \V\ is maximal (see Appendix C of iKrainovic et aU 
l|2006ll l. Triaxial systems can have complicated velocity 
fields, because of the possible superposition of orbits with 
angu lar momenta with respect to the minor and the major 
axis (|van den Bosch et al.ll2008l ). In a perfect axisymmetric 
galaxy there can be no kinematic misalignment. Indeed all 
fast rotators in the SAURON sample have misalignments 

< 5°, whereas slow rotators can have any kinematic mis- 
alignment up to 90°. Our example remnants are also quite 
diverse in their behaviour of kinematic misalignment with 
projection angle (Fig[5]). The oblate remnant as it is dom- 
inated by minor axis tubes has very few projections where 

7^ 0. The triaxial remnant has a much broader distribu- 
tion, which peaks at low misalignments, but the maps of 
this remnant show little rotation, such that the misalign- 
ment has larger error bars and is not so meaningful. The 
prolate remnant, however, shows clearly a strong misalign- 
ment, i.e. > 20, for most of the projections, as both tube 
orbit classes revolving around the major as well as the minor 
axis, are present in this remnant. 



5 SLOW AND FAST ROTATORS 

The determination of the A-parameter is important in sofar 
as ellipticals can be divided into two broad sub-classes, so 
called slow and fast rotators (EMQ7). The (empirical) cut 
which separates these classes has been made at Xr = 0.1. 
These seem to be genuinely intrinsically different galaxy 
types with little overlap through projection effects. Many 
slow rotators exhibit e.g. large, old and fast rotating kine- 
matically decoupled components (KDCs), which however 
are not affecting too much the overall angular momentum 
balance, as the stars which move on high angular momen- 
tum orbits are not located in the center, such that the Xr 
as determined from the map is rather low. In general maps 
of slow rotators resemble more the kinematic maps derived 
from equal-mass merger remnants. The formation of KDCs 
requires the inclusion of a gaseous component in the simu- 
lation. Fast rotators are on the other hand more likely to be 




Figure 6. Distribution of for coUisionless merger remnants of 
different merger mass ratios. The region where kinematic maps 
are classified as slow rotators is shaded. Only binary mergers of 
equal-mass progenitors seem to produce slow rotators in sizeable 
fractions. 

un-equal mass mergers. For a detailed discussion on the ori- 
gin of kinematical features in 2-dimensional LOSVD maps 
we refer the reader to J07. 



5.1 The Role of Merger Mass Ratios 

The most fundamental parameter which influences the out- 
come of a merging event is the amplitude of the fluctua- 
tions of the gravitational potential during the merger. How 
strongly the phase space is rearranged after a merger will be 
mainly determined by the ratio of the masses of the merging 
objects (,Naab fc Burkert 2003). An equal-mass merger will 
be more effective in redistributing the energies of the stars 
in the progenitor galaxies than an un-equal mass merger. It 
is true that secondary factors like infall velocity and merging 
geometry have also an impact on the final shape of the rem- 
nant, but the mass ratio of the merger is the dominant fac- 
tor. We want to quantify when a given projection of a given 
remnant is classified as a slow or a fast rotator and we use 
coUisionless disc-disc merger remna nts of mergers wi t h mas s 
ratios of 1:1, 2:1, 3:1 and 4:1 from iNaab fc BurkertI (|2003l ). 
In Fig[6]the Xr distribution is shown for each of these merger 
subsamples. We see, as expected, that the fraction of slow 
rotators for a given merger mass ratio is decreasing with in- 
creasing mass ratio. Xr peaks at 0, 0.1, 0.27 and 0.4 for the 
mass ratios 1:1 to 4:1. If one would assume that binary disc- 
merging is the main formation channel of elliptical galaxies 
then only equal or near equal-mass merger could produce 
them. However, this picture is probably too simple as ellip- 
tical galaxies can have more complicated merger histories i n 
a cosmological context (|Meza et al.ll2003HNaab et al.ll2007l '). 
Multiple minor mergers, e.g. three 3:1 mergers might also 
be a viable way to form slowly rotating ellipticals, which 
resemble closely the outcome of a single equal-mass merger, 
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Figure 4. Left: Case study of projectional relation between Xn and photometric properties, i.e. ellipticity and isophotal shape parameter 
for an oblate (red), a prolate (black) and a triaxial remnant (blue). The triaxial remnant has a Ri not varying with projection. 
Ellipticity and Xr are good proxies for inclination angle for the oblate remnant, which is nearly axisymmetric. The prolate remnant can 
have more than one value for a given ellipticity owing to a complicated orbital structure. Right: Histogram of ellipticities and A^ of 
the projectional study. The oblate remnant is a fast rotator for almost all viewing angles, while the triaxial remnant is a slow rotator for 
all viewing angles. The prolate remnant is equally probable to be observed as a fast or slow rotator. 



Merger Mass Ratio 


Slow Rotator Fraction 


1:1 


0.75 


2:1 


0.10 


3:1 


0.03 


4:1 


0.02 



Table 1. Dependence of the slow rotator fraction on the merger 
mass ratio. Only equal-mass merger are very effective in removing 
angular momentum in a single merging event. 

as was shown recently bv lBournaud et all ()2007h . All the re- 
sults (except for Sec l5.4p we present in this paper are valid 
for one generation of merging only. The examination of the 
impact of complex hierarchical merging histories is beyond 
the scope of this paper. 

5.2 Confusion Rate 

ICappellari et al.] (|2007l ) showed that slow and fast rotators 
are indeed intrinsic fast and slow rotators and not randomly 



projected fast rotating galaxies which happen to appear non- 
rotating. We address this question in Fig. [7] where we exam- 
ine the fraction of projections for each of the collisionless 
remnants which have Xr < 0.1. There is an isolated sample 
of remnants which have a probability of greater than 90% to 
be classified as slow rotators. The remaining remnants have 
a probability of less than 40% to be identified as slow rota- 
tors. The mergers which are located at the extreme ends of 
the distribution are the unambiguous cases they are (almost) 
always either fast rotators or slow rotators. There is how- 
ever a small population of remnants which have a reasonable 
probability to be identified as slow rotators although they 
have a non-zero angular momentum content (shaded zone in 
FiglZl left) . These are on average the most prolate of all rem- 
nants (same Figure, right plot). The confusion arises from 
the presence of major and minor axis tubes in the remnants, 
i.e. there are two face-on projections under which the angu- 
lar momentum content is perpendicular to the line-of-sight 
and cannot be observed. But there also remnants which are 
mildly to strongly triaxial (0.3 < T < 0.5), which are hot but 
regularly rotating systems, i.e. with low kinematic misalign- 
ment. Almost all of these systems resulted from 2:1 merg- 
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Figure 5. Left: Case study of kinematical misalignment ^ of the same remnants as in Fig|4] Relation between Aj; and kinematic 
misalignment. Right: Distribution of kinematical misalignments from 200 random projections. The OBLATE remnant has very few 
projections with ^ > 5° , while the PROLATE remnant distribution of misalignment angles peaks at 40° . The TRIAXIAL remnant has 
a less broad distribution of 'I' values. 



ers. We can calculate now the global probability to falsely 
classify a rotating galaxy as a slow rotator in our sample 
and find Pconfuso = 0.046. This is probably an upper limit 
as prolate remnants are overabundant in our sample, and 
are suppressed if a significant am ount of gas is included in 
the simulations (iNaab et al.ll2006i ). This is consistent with 
the SAURON sample where not a single galaxy has been 
observed with rotational patterns as, e.g. our PROLATE 
model, which we presented in Sec|3] We further note that 
the maximal Aj? is a good tracer for the triaxiality of our 
remnants (Fig[7lright), albeit with an increasing spread to- 
wards lower Aij. Particular spin alignments of the progenitor 
discs can lead to interesting outliers from this correlation, 
e.g. we find a slow rotator with T ~ 0. This merger remnant 
originated from a planar merger of discs with anti-aligned 
spins. The counter-rotation of both disc components cause a 
net zero rotation. While this is a possible formation scenario, 
most slow rotators are triaxial box dominated remnants. 



5.3 The Influence of Dissipation 

Late-type spiral galaxies are normally not purely collision- 
less systems but have a sizable fraction of gas. We there- 
fore compare two sets of mergers: disc-disc mergers with gas 
and star formation, and coUisionless disc-disc mergers. Each 
set consists of 16 1:1 and 32 3:1 mergers, which formed on 
identical merging geometries. As mentioned before we take 
200 observations of each remnant at random viewing angles, 
which results in 9600 maps for each set of simulations. We 
determine for each set the fraction of slow and fast rotators, 
as well as standard photometric parameters such as elliptic- 
ity and the isophotal shape parameter a^. In Fig.|S]we show 
the distributions of ellipticities and 04 for slow and fast rota- 
tors from mergers which formed with and without gas. Maps 
classified as slow and fast rotators show distinct photomet- 
ric properties. Slow rotators in general have smaller elliptic- 
ities and have smaller 04 values while fast rotators are more 



elliptical and discy. It is apparent that the fast rotator pop- 
ulation hardly changes its photometric properties when gas 
is included in the simulation. The impact on slow rotators, 
however, is more visible as the ellipticity distribution is now 
more ske wed towards e = 0, while the 04 d istribution peaks 
now at (|Cox et al.ll2006llNaab et al.ll2006l '). The mean slow 
rotator properties are summarized in Table [2j which also 
confirm that most slow rotators originate from equal-mass 
mergers. The results from the SAURON survey showed that 
slow and fast-rotating ellipticals are located at certain posi- 
tions in the A — e and in the A — a4 plane. For a better com- 
parison we also show the distribution of merger remnants 
in those parameter spaces (Fig.|9]). In general the parameter 
space bracketed by the SAURON observations is reproduced 
well by the merger remnants. There are, however, some dif- 
ferences. For example we see that the coUisionless mergers 
have too high e for a given Xr which is alleviated, especially 
for slow rotators, by the inclusion of gas. We also indicated 
the progenitor galaxy in Fig|9] We already showed in Sec l5.ll 
that the amount of rotational support always decreases af- 
ter a merger event. Very high merger mass ratios will lead 
to merger remnants which would lie closer to the progenitor 
limit, but we have not performed 10:1 or even higher mass 
ratios. Therefore we do not find projections of merger rem- 
nants with Xr > 0.6. It is remarkable though that some of 
the SAURON fast rotators have (A_R,e)-values very close to 
the progenitor galaxy and are, at least in this sense, hardly 
distinguishable from early-type spiral galaxies. 

The situation is different in the 04 — A-plane where the 
merger remnants are only partly located in the parameter 
space occupied by the SAURON galaxies. CoUisionless slow 
rotators are too boxy while fast rotators are on average too 
discy compared to SAURON galaxies. While gas removes the 
boxy projections of slow rotators it does not help to create, 
as expected, boxy projections in fast rotators. Probably the 
origin of boxiness in the fast rotating SAURON galaxies is 
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Figure 7. Left: Distribution of mergers according to tlieir probability of being classified as a slow rotator (Aj; < 0.1). Right: Correlation 
between the maximal and the intrinsic triaxiality T of the merger remnants. Open symbols indicate merger remnants which belong 
to the shaded area from the plot on the left side, i.e. have a significant number of misclassified projections. 
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Table 2. Slow Rotator Properties. Main results of the classifica- 
tion according to the A^-parameter. Every projection which has 
a A < 0.1 is classified as a slow rotator. Dissipation mainly in- 
fluences the ellipticity, while the slow rotator fraction is almost 
unaffected. 
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Figure 8. Distribution of ellipticitios and isophotal shape pa- 
rameter for the slow rotator (left plots) and the fast rotator (right 
plots) population. The collisionless 1:1 and 3:1 remnants are indi- 
cated by black solid lines. Fast and slow rotators have clearly sep- 
arate photometric properties. Fast rotators have high ellipticities 
and discy isophotes. Histograms for dissipatively formed merger 
remnants are overplotted in (dashed) red. Fast rotators are very 
little affected while slow rotators with gas have significantly lower 
ellipticities, i.e. are rounder. 



diflterent from the source of boxiness in our remnants, e.g. 
bars. 

As we discussed earlier, kinematical misalignment is 
also a distinguishing feature between fast and slow rotators. 
Again we divide up the simulation sample into fast and slow 



rotators and measure the angle ^' between kinematic and 
photometric position angles. The observational trends are 
reproduced remarkably well: fast rotators are very unlikely 
to show kinematic misalignment, while slow rotators show 
a flat distribution of misalignment angles, in other words 
there is no preferred alignment angle for slow rotators (see 
FigHnj- The low misalignment angles for fast rotators are 
not very surprising, because we already start out with a per- 
fectly aligned system. Unequal mass merger simply do not 
destroy the initial alignment. Indeed observed fast rotators 
are also remarkab ly aligned. Only 2 out of 50 fast rotators 
have * > 5° (see ICappellari et al.|[2007l . Fig.l2). The cru- 
cial question to explain the origin of fast rotators is rather 
how the axisymmetric progenitor systems are produced in 
the first place, rather than the subsequent merging. 

Gas changes the picture only slightly. There are some- 
what fewer projections with high misalignment angles for 
fast rotators and slow rotators are slightly better aligned. 
In general the inclusion of only 10% of gas in the merging 
process is not necessary to explain the kinematic misalign- 
ment of slow and fast rotators. However, the inclusion of 
high fractions of gas would probably change this picture. 



10 R. Jesseit, M. Cappellari, T. Naab, E. Emsellem and A. Burkert 




10. 

17. 

|33. 
50. 

^ i67. 

83, 
100. 



0.0 0.2 0.4 0.6 0.8 



1.0 
OA 



Collisionless Mei'gei's 
★ SAURON Galaxies 
Progenitor Galaxy 



1.0 

0.8 




Mer'gers with Gas + SF 




Figure 9. Left Column: Two-dimensional probability distribution of coUisionless 1:1 and 3:1 remnants in — 04, respectively in the 
Afl — e plane. The SAURON galaxies are indicated by black stars. High A^ galaxies are resembling more the progenitor galaxy than 
3:1 merger remnants. Boxy fast rotators are not formed in mergers, but are probably barred galaxies. Right Column: The same plots 
but this time with merger remnants which formed with gas. The photometric properties are stronger affected than Ajj. Slow rotating 
remnants agree well with SAURON galaxies. 



5.4 Dry Mergers 

The most massive galaxies in the universe are probably 
not form ed from a single genera tion of binary spiral-spiral 
mergers (iNaab fc Ostriker l2007h . because there are simply 
not enough massive late-type galaxies to account for the 
stellar mass of galaxies with luminosities, e.g., of 4 L* . 
Their merging history included probably more than one 
major merger and semi-analytic mode ling showed that the 
last major merger was p robably dry (jKhochfar fc BurkertI 
l2003llHopk ins et al] 120071 ). We have a small sample of re- 
mergers of coUisionless merger remnants to test the assump- 
tion if slow rotators are the end product of dry merging. 
We analyzed 6 equal-mass mergers of discy remnants and 6 
equal-mass mergers of boxy remnants. The majority of the 
maps (about 68 %) of the dry mergers are indeed slow rota- 



tors. As we have a very limited sample of merger remnants, 
we compare them only with the slow rotators which orig- 
inated from 1:1 coUisionless disc mergers, which is admit- 
tedly the overwhelming majority. In Fig. [11] we show that 
the ellipticity distribution from elliptical-elliptical mergers 
is hardly distinguishable from equal-mass disc-disc mergers. 
However, they are on average more boxy than r emnants of 
disc mergers in agreement with lNaab et al.l (j2006l ). The slow 
rotators in the SAURON sample are indeed hardly boxy. As 
mentioned before, the SAURON sample is not representa- 
tive of the shape distribution of all elliptical galaxies due 
to the selection procedure of galaxies. There is especially a 
dearth of triaxial, boxy galaxies which makes it difficult to 
exclude or confirm the dry merging ori gin of massive ellipti - 
cals with the present sample (see also iBurkert et al.l (|2007f ) 
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Figure 10. Histograms of kinematic misalignment angle ^ mea- 
sured from 2-dimensional kinematical maps of the merger rem- 
nants. The overwhelming majority of maps identified as fast ro- 
tators have no kinematic misalignment, while slow rotators have 
no preferred misalignment angle. Dissipation changes the results 
only slightly, i.e. the maps slow rotators have appear slightly more 
aligned. 
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Figure 11. Distribution of photometric properties for slow ro- 
tators from coUisionless 1:1 elliptical and 1:1 disc mergers. The 
ellipticity distribution of both samples are similar while elliptical 
merger remnants are more boxy. 



for a more detailed discussion of the origin of the slowly 
rotating SAURON ellipticals). 



6 DISCUSSION AND CONCLUSIONS 

We calculated the A_R-parameter from line-of-sight velocity 
maps of 48 remnants of coUisionless disc mergers, 48 rem- 
nants of disc mergers with a dissipative component and 12 
re-mergers of coUisionless remnants, according to the proce- 
dure laid out in EM07. Every remnant is projected 200 times 
randomly and then treated as a mock galaxy. We first inves- 
tigated how the intrinsic properties of the merger remnants 
manifest itself in the maps and how much can be deduced of 
the intrinsic structure by calculating Xr. We then identified 
a mock galaxy as a fast or slow rotator and compare their 
photometric and kinematic properties like 04, e and kine- 
matic misalignment <& to the properties of SAURON galax- 
ies. We examined the influence of gas and merging mass 
ratio on the statistical properties of slow and fast rotators. 
We report the following findings: 

• The line-of-sight angular momentum parameter Xr is a 
rather robust indicator of the intrinsic angular momentum 
content of a galaxy. Xr stays close to its maximum value for 



a wide range of viewing angles as it captures the balance 
between ordered (line-of-sigh velocity) and unordered mo- 
tion (line-of-sight velocity dispersion) well. In principle the 
projected angular momentum content could also be deduced 
from the two-dimensional velocity field alone, but informa- 
tion over the entire field of the galaxy would be needed, 
which is normally not available. We conducted therefore 
tests to explore the deviation from the intrinsic angular mo- 
mentum, observing galaxies with apertures of limited size. 

• The projectional behaviour of Xr is determined by the 
intrinsic orbital structure. Most fast rotators consist of disc- 
like minor axis tubes and have A_r > 0.1 for almost all pro- 
jections. Likewise slow rotators consist of box orbits and 
have A_R < 0.1 for all projections. Prolate remnants have 
a more complicated projectional behaviour as they consist 
of major and minor axis tubes. There is however not a sin- 
gle galaxy in the SAURON survey with significant rotation 
around the major axis and we slightly overproduce these sys- 
tems in our merger sample. Even so we determine the mis- 
classification probability for the whole coUisionless merger 
sample to only 4.6 %. We therefore confirm the robustness 
of the slow/fast rotator classification scheme for binary disc 
merger remnants. 

• The spread of Xr is naturally reproduced in the binary 
merger scenario. However their photometric properties have 
some important differences. Slow rotators from coUisionless 
mergers have too high ellipticities. The problem is alleviated 
by the inclusion of a modest fraction of gas, which does not 
modify the slow rotator fraction but produces remnants with 
lower ellipticities. Therefore gas, maybe counter-intuitively, 
can be important for the formation of slow rotators. 

• The destruction of ordered rotation is only effective in 
equal-mass mergers. The violence of the merger seems to 
be the decisive factor, i.e. with increasing mass-ratio the 
percentage of formed slow rotators drops very rapidly. This 
is true at least for the binary merger scenario, very large gas 
fractions and multiple mergers will complicate this picture. 

• We tested the amount of kinematic misalignment for 
slow and fast rotators separately. Again the agreement with 
observations is very good. Slow rotators can have any mis- 
alignment with equal probability, while fast rotators show 
kinematical alignment for almost all projections. The pic- 
ture changes with the inclusion of gas only insofar as the 
strongest misalignments vanish for slow rotators. Higher gas 
fraction will probably lead to a further decrease in observed 
misalignments. 

• Slow rotators in the SAURON sample are not strongly 
boxy and rather round, therefore the dry mergers are a poor 
fit to them. This is somewhat surprising, as the projection 
and mass ratio independence o f boxiness was seen as an asset 
of this formation mechanism l|Naab et al.ll2006l ). To assess 
the role of dry mergers adequately probably a larger sample 
of eUipticals, like ATLAS^°, is needed. 

• Fast rotators with boxy isophotes are absent from our 
sample. Some early-type SOs in the SAURON are barred 
galaxies, and such strong bars seem not to be formed in our 
particular merger sample. Our sample of course has its limits 
and it is possible that the bulge to disc ratio of our progeni- 
tor galaxies is too high or the initial bulges are not rotating 
or probably both, to compare to this particular property of 
the SAURON galaxies. 

• The fast rotators with the highest Xr values in the 
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SAURON sample resemble our progenitor galaxy more than 
a merger remnant, which leads to the question, if major 
merging played a significant role in their formation at all. 
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